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The framework of relativistic energy density functionals is extended to include correlations related 
to the restoration of broken symmetries and to fluctuations of collective variables. The generator 
coordinate method is used to perform configuration mixing of angular-momentum projected wave 
functions, generated by constrained self-consistent relativistic mean-field calculations for triaxial 
shapes. The effects of triaxial deformation and of if-mixing is illustrated in a study of spectroscopic 
properties of low-spin states in ^*Mg. 
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I. INTRODUCTION 

Among the microscopic approaches to the nuclear 
many-body problem, the framework of nuclear energy 
density functionals (EDF) is the only one that can 
presently be used over the whole nuclear chart, from rel- 
atively light systems to superheavy nuclei, and from the 
valley of /3-stability to the particle drip- lines [ll-Q . Mod- 
ern energy density functionals provide the most complete 
and accurate description of structure phenomena related 
to the evolution of shell structure in medium-mass and 
heavy nuclei, e.g. the appearance of new regions of de- 
formed nuclei, shape coexistence and shape transitions. 

In practical implementations the EDF framework is 
realized on two specific levels. The simplest implemen- 
tation is in terms of self-consistent mean-field models, 
in which an EDF is constructed as a functional of one- 
body nucleon density matrices that correspond to a sin- 
gle product state - Slater determinant of single-particle 
or single-quasiparticle states. This framework can thus 
also be referred to as single reference (SR) EDF. In the 
self-consistent mean-field approach the many-body prob- 
lem is effectively mapped onto a one-body problem, and 
the exact EDF is approximated by a functional of pow- 
ers and gradients of ground-state nucleon densities and 
currents, representing distributions of matter, spins, mo- 
mentum and kinetic energy. In principle the SR nuclear 
EDF can incorporate short-range correlations related to 
the repulsive core of the inter-nucleon interaction, and 
long-range correlations mediated by nuclear resonance 
modes. The static SR EDF is characterized by symme- 
try breaking - translational, rotational, particle number, 
and can only provide an approximate description of bulk 
ground-state properties. To calculate excitation spectra 
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and electromagnetic transition rates in individual nuclei, 
it is necessary to extend the SR EDF framework to in- 
clude collective correlations related to the restoration of 
broken symmetries and to fluctuations of collective co- 
ordinates. Collective correlations are sensitive to shell 
effects, display pronounced variations with particle num- 
ber, and cannot be incorporated in a SR EDF. On the 
second level that takes into account collective correlations 
through the restoration of broken symmetries and con- 
figuration mixing of symmetry-breaking product states, 
the many-body energy takes the form of a functional of 
all transition density matrices that can be constructed 
from the chosen set of product states. This level of im- 
plementation is also referred to as multireference (MR) 
EDF framework. 

In recent years several accurate and efficient models 
and algorithms have been developed that perform the 
restoration of symmetries broken by the static nuclear 
mean field, and take into account fluctuations around 
the mean-fleld minimum. The most effective approach 
to configuration mixing calculations is the generator co- 
ordinate method (GCM) [1, With the simplifying as- 
sumption of axial symmetry, GCM configuration mixing 
of angular-momentum, and even particle-number pro- 
jected quadrupole-dcformcd mean-field states, has be- 
come a standard tool in nuclear structure studies with 
Skyrme energy density functionals [l|, the density- 
dependent Gogny force 0], and relativistic density func- 
tionals [1, 0. A variety of structure phenomena have 
been analyzed using this approach. For instance, the 
structure of low-spin deformed and superdeformed col- 
lective states fn)l -[l2|. shape coexistence in Kr and Pb 
isotopes [l^ Il4{ . shell closures in the neutron-rich Ca, 
Ti and Cr isotopes [l5| and shape transition in Nd iso- 
topes [il,[i3. 

Much more involved and technically difficult is the de- 
scription of intrinsic quadrupole modes including triaxial 
deformations. Intrinsic triaxial shapes are essential for 
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the interpretation of interesting collective modes, such 
as chiral rotations [H, [l^ and wobbling motion poj . 
The inclusion of triaxial shapes can dramatically reduce 
barriers separating prolate and oblate minima, leading 
to structures that are soft or unstable to triaxial dis- 
tortions [llj. Such a softness towards dynamical 7- 
distortions will give rise to the breakdown of the de- 
selection rule in electromagnetic transitions of high-spin 
isomers [l^. It may also has important influence on 
the electric monopole transition strength B{EO : — > 

0+) p. 

Only very recently a fully microscopic three- 
dimensional GCM model has been introduced p^ ]. 
based on Skyrme mean-field states generated by triax- 
ial quadrupole constraints that are projected on parti- 
cle number and angular momentum and mixed by the 
generator coordinate method. This method is actually 
equivalent to a seven-dimensional GCM calculation, mix- 
ing all five degrees of freedom of the quadrupole op- 
erator and the gauge angles for protons and neutrons. 
In this work we develop a model for configuration mix- 
ing of angular-momentum projected triaxial relativistic 
mean-field wave functions. In the first part, reported in 
Ref. [25| . we have already considered three-dimensional 
angular-momentum projection (3DAMP) of relativistic 
mean-field wave functions, generated by constrained self- 
consistent mean-field calculations for triaxial quadrupole 
shapes. These calculations were based on the relativistic 
density functional PC-Fl [2^, and pairing correlations 
were taken into account using the standard BCS method 
with both monopole and zero-range d interactions. Cor- 
relations related to the restoration of rotational symme- 
try broken by the static nuclear mean field, were ana- 
lyzed for several Mg isotopes. Here we extend the model 
of Ref. [2^ , and perform GCM configuration mixing of 
3DAMP relativistic mean-field wave functions. 

In Section |lT] we introduce the model, briefly outline 
the relativistic point-coupling model which will be used 
to generate mean-field wave functions, and describe in 
detail the procedure of configuration mixing of angu- 
lar momentum projected wave functions. In Section IIIII 
the 3DAMP-I-GCM model is tested in illustrative calcu- 
lations of the low-energy excitation spectrum of ^^Mg. 
Section HVl summarizes the results of the present investi- 
gation and ends with an outlook for future studies. 



II. THE 3DAMP+GCM MODEL 

The generator coordinate method (GCM) is based on 
the assumption that, starting from a set of mean-field 
states I $((/)) which depend on a collective coordinate 
q, one can build approximate eigenstates of the nuclear 
Hamiltonian 



Detailed reviews of the GCM can be found in 
In the present study the basis states \^{q)) are Slater 
determinants of single- nucleon states generated by self- 
consistent solutions of constrained relativistic mean-field 
(RMF) -I- BCS equations. To be able to compare the- 
oretical predictions with data, it is of course necessary 
to construct states with good angular momentum. Thus 
the trial angular-momentum projected GCM collective 
wave function \ an eigenfunction of and Jz, with 

eigenvalues J{J + l)fi^ and Mfi, respectively, reads 
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\JMK+,q) (2) 



where a — 1,2, ■ ■ ■ labels collective eigenstates for a given 
angular momentum J. The details of the 3D angular- 
momentum projection are given in Ref. [25| . here we 
only outline the basic features. Because of the D2 and 
time-reversal symmetry of a triaxially deformed even- 
even nucleus, the projection of the angular momentum 
J along the intrinsic z-axis {K in Eq. ([5]) ) takes only 
non-negative even values: 



0, 2, • • • , J for J mod 2 
2,4, • • • , J - 1 for J mod 2 = 1 



(3) 



The basis states \JMK+,q) are projected from the in- 
trinsic wave functions \^{q)): 

\JMK+, q) = [P-U + {-lYPii^Kmq)), (4) 

where P'lij^ is the angular-momentum projection opera- 
tor: 



2 / + 1 

PiiK = I dnDi*j,in)R{n) 



87r2 



(5) 



H. denotes the set of three Euler angles: {(f>,9,ip), and 
dfl ~ d(j) sin 9 d9 dtp. Dfjj^{n) is the Wigner D-function, 
with the rotational operator chosen in the notation of 
Edmonds [13 : R{^) = e">''^- e^^^y e'"^^- . The set of in- 
trinsic wave functions \<^{q)), with the generic notation 
for quadrupole deformation parameters q = (/3,7), is 
generated by imposing constraints on the axial (720 and 
triaxial 922 mass quadrupole moments in self-consistent 
RMF+BCS calculations. These moments are related to 
the Hill- Wheeler [H] coordinates ;3 (^ > 0) and 7 by the 
following relations: 



(1) 
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-|- (2z2 -x^-y^)^ ^ARl(3 COS 7, (6a) 
loTT 47r 



^y^y')-l^Pl^,P^^^^ (6b) 
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where Rq = 1.2A^^^ fm. The total mass quadrupole mo- 
ment (7,„ reads: 



IQtt 



(7) 



The calculation of single-nucleon wave functions, ener- 
gies and occupation factors starts with the choice of the 
energy density functional (EDF) . As in our previous anal- 
ysis on collective correlations in axially deformed nuclei 
@, Q , and in the first part of this work [25| , the present 
illustrative calculation is based on the relativistic func- 
tional PC-Fl (point-coupling Lagrangian) [26j : 



Ebmf ^ y £viMFir) (8) 

= ^ dr vl ipkir) (-17 V + m) tpkir) 
k 

+ -^.iTvi^Tv)t. + —jrv^ijTv),, 

+ -I^Pts +—Pts^Pts + e — —PvA 

where ipki^) denotes a Dirac spinor. The local isoscalar 
and isovector densities and currents 



Psir) 
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fe>0 




PTs{r) 


= Xl^fc V'fc(r-)T3V'fc(r) , 


(9b) 




fc>0 




f{r) 
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= ^kir)-/^T3^Pk{r) , 


(9d) 



fc>0 

are calculated in the no-sea approximation, i.e., the sum- 
mation runs over all occupied states in the Fermi sea. 
The occupation factors u| of each orbit are determined 
in the simple BCS approximation, using a (5-pairing force. 
The pairing contribution to the total energy is given by 

EpM^^,K*] = - E T / d'r^*Ar)nr{r). (10) 

r—n.p 

where Vt is a constant pairing strength, and the pairing 
tensor K{r) reads 



K{r) = -2^ /fcUfeDfc|-0fc(r)|^ 



(11) 
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factors fk , determined by a Fermi function in the single- 
particle energies e^: 



fk 



1 



1 -I- exp[(efc -ep - AEr)/pT] 



(12) 



ef is the chemical potential determined by the constraint 

on average particle number: ($((7)|iVr|$((j')) = Nr- The 
cut-off parameters AEt and Pr = AEr/10 are chosen in 

such a way that 2 = N^. + 1.65Nr^'^, where is 

the number of neutrons (protons) j29| . 

The weight functions fi^{q) in the collective wave 
function Eq. ([2]) are determined from the variation: 



JM\ \^JM 







(13) 



i.e., by requiring that the expectation value of the energy 
is stationary with respect to an arbitrary variation Sf^^. 
This leads to the Hill- Wheeler- Griffin (HWG) integral 
equation: 



A">0 

, (14) 

where ^ and ^ are the angular-momentum projected 

GCM kernel matrices of the Hamiltonian and the norm, 
respectively. With the generic notation ^ = or J^, 
the expression for the kernel reads: 

+ {-iyOi_j,,{q,q') + i-iyOiKK,iq,q% 

(15) 



where for the operator O = 1 or iJ: 



Oij,,{q,q')^{<i>{q)\dP^^,mq')}, 



(16) 



and Akk' = 1/[(1 + Sko){1 + Sk'o)]. 

The overlap {^{q)\HR\^{q')) can be evaluated in co- 
ordinate space, and we rewrite the hamiltonian kernel 
Hj^j^,{q, q') in the following form: 

H^K' (9, q') = J drH^KK' (r; q, q') , (17) 



where 



2 7+1 

HiiK'{r:q,q')^^-^ / d^lDi\,n{r-q,q'-^)n{q,q'-^) 



The pairing window is constrained with smooth cutoff 



The norm overlap n{q, q'; il) is defined by: 
n{q,q';n) = mq)\R{nmq')) 



(18) 



(19) 
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The calculation of the overlap matrix elements state-independent value of the Lagrange parameter, that 



7i{r;q,q';il) requires the explicit form of H. So far we 
have implicitly assumed that the system is described by 
a Hamiltonian. However, for energy density function- 
al this is strictly valid only if the density dependence 
can be expressed as a polynomial of p. By using prod- 
uct wave functions, a density functional can formally be 
derived from a Hamiltonian that contains many-body in- 
teractions. A prescription based on the generalized Wick 
theorem j30l | states that the Hamilton overlap matrix el- 
ements have the same form as the mean field functional, 
with the intrinsic single particle density matrix elements 
replaced by the corresponding transition density matrix 



elements |3l| . In this work we employ the relativistic 
point-coupling model PC-Fl [2^, which contains powers 
of the scalar density ps up to fourth order, and therefore 
the above prescription can be applied. For a detailed dis- 
cussion of open problems we refer the reader to Ref . [12] , 
and references cited therein. 

Consequently, 7^(r; g'; fi) has the same form as the 
mean- field functional £rmf{i') in Eq. ^ provided the 
intrinsic densities and currents are replaced by transi- 
tion densities and currents. Further details about the 
calculation of the norm overlap n{q, q'; fl) and transition 
EDF H(r; q, q'; il) can be found in Ref. [25|. 

The basis states |$(q)) are not eigenstates of the pro- 
ton and neutron number operators Z and N. The adjust- 
ment of the Fermi energies in a BCS calculation ensures 
only that the average value of the nucleon number oper- 
ators corresponds to the actual number of nucleons. It 
follows that the wave functions are generally not 

eigenstates of the nucleon number operators and, more- 
over, the average values of the nucleon number operators 
are not necessarily equal to the number of nucleons in 
a given nucleus. This happens because the binding en- 
ergy increases with the average number of nucleons and, 
therefore, an unconstrained variation of the weight func- 
tions in a GCM calculation will generate a ground state 
with the average number of protons and neutrons larger 
than the actual values in a given nucleus. In order to 
restore the correct mean values of the nucleon numbers, 
we follow the standard prescription [33;, i3^] , and modify 
the HWG equation by replacing 7i{r; q, q'; ft) with 

U'{r;q,q';n) - H{r; q, q' ;n) - \[Z (r; q, q' ■,n) - Zo] 
-Xn[N{r;q,q';n)-No], (20) 

where Zq and A'o are the desired proton and neutron 
numbers, respectively. Z{r; q, q'; 51) and N{r; q, q'; fi) are 
the transition vector densities in r-space for protons and 
neutrons, respectively. The Lagrange parameters At=p,ti 
are in principle determined in such a way that each AMP 
GCM collective state has the correct average particle 
number. In that case, however, the Lagrange parameters 
At- will be state dependent and, as a consequence, the 
orthonormality of the states |4';^*'^)s is no longer guaran- 
teed. In Ref. |34| a simple ansatz was introduced for a 



is the value of Ar=p,n was chosen to be the mean BCS 
Fermi energy, determined by averaging over the collective 
variable q. The average particle numbers in the resulting 
AMP GCM states differ only slightly from the desired 
correct values. In the present model we take the same Ar 
values as those in the mean-field calculation, i.e. Xt{q) 
for the diagonal terms {q' = q), and [A,-(g) -I- At-('z')]/2 for 
the off-diagonal ones {q' 7^ q) in T-L'{r;q,q';n). We find 
that with this prescription the average particle numbers 
for low-lying excitation states are in excellent agreement 
with those obtained by taking the A7- value averaged over 
the collective variable q. 

The domain of quadrupole deformation parameters 
q = {f3, 7) is discretized, and the HWG integral equation 
is transformed into a matrix eigenvalue equation. The 
corresponding kernels ^KK'il^l') have to be calculated 
between all pairs of mesh points in q space. In the cur- 
rent version of the model the full space K ^q\s a. direct 
product of the if-subspace and the g-subspace, with di- 
mension D — {J + 2)nq/2 for even J or D = (J — l)ng/2 
for odd J. riq is the number of points on the mesh in q- 
space, and J the total angular momentum. Correspond- 
ingly, the kernels ff^K' il' 1') 
number K and the value of (/ 

space K ^q can be determined as shown in Fig. [TJ 

The first step in the solution of the HWG matrix eigen- 
value equation is the diagonalization of the norm overlap 
kernel ,yV'^{i,j) 



The quantum 
7) at each point of the full 



(21) 



Since the basis functions \^{q)) are not linearly inde- 
pendent, many of the eigenvalues are very close to 
zero. They correspond to "high momentum" collective 
components, i.e., the corresponding eigenfunctions u^l{i) 
arc rapidly oscillating in the g-space but carry very lit- 
tle physical information. However, due to numerical un- 
certainties, their contribution to the matrix elements of 
the collective Hamiltonian can be large, and these 
states should be removed from the basis. Therefore, a 
small positive constant C is introduced so that states 
with n^/n^^^ < ^ are excluded from the GCM basis, 
where n^^^, is the largest eigenvalue of the norm kernel. 
From the remaining states, also called ^^natural states^', 
one builds the collective Hamiltonian 



njj 
Ttuj 



1 



1 



ui{{).^\i,3)ui{j) , (22) 



which is subsequently diagonalized 



(23) 



The solution of Eq. ([23| determines both the energies 
and the amplitudes f^^ {q) of collective states with good 
angular momentum l^*^^^) 
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J: Even 
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FIG. 1: (Color online) Distributions of the K quantum number (K) and q values (iq) in the full K space. The dimension 
of the direct product of the -RT-subspace and the g-subspace is _D = ( J + 2)nq/2 for even J, or D — {J — l)nq/2 for odd J. Uq 
is the number of points on the mesh in g-space, and J is the total angular momentum. 



(24) 



The weight functions f^^ (q) are not orthogonal and can- 
not be interpreted as collective wave functions for the de- 
formation variables. The collective wave functions 5^(j) 
arc calculated from the norm overlap eigcnstates: 



(25) 



g^{i) are orthonormal and, therefore, |(7Q(j)P can be in- 
terpreted clS &- probability amplitude. 

The center-of-mass (cm.) correction is defined by: 



p2 

J A/ I cm |yj;.7Af^ 



'2mA 



i] KK' 



2mA 



(26) 



The projected overlap matrix elements 
($((7i)|P^„jPj^^,|$(gj)) are treated in zeroth order 
of the Kamlah approximation, i.e. considering the fact 
that {^{qi)\^{qj)) is sharply peaked at qi = qj, the 
projected matrix elements are approximated by the 
unprojected ones and 



1 



2mA 



(27) 

were Ecm{q) is the cm. correction evaluated for the in- 
trinsic wave functions \'^{q)), 

1 



2mA 



($(q)|PLl$(9)) , 



(28) 



where m is the nucleon mass, and A is the number of 
nucleons. Pcm = X]f total momentum. The 

energy of the collective state j^E';^) is, therefore, given by 



(29) 



Once the amplitudes fa^ (q) of nuclear collective wave 
functions l^'^^^) are known, it is straightforward to cal- 
culate all physical observables, such as the electromag- 
netic transition probability, spectroscopic quadrupole 
moments and the average particle number. The B(E2) 
probability for a transition from an initial state {Ji,ai) 
to a final state (J/, a/) is defined by 



{Jf^1f\\Q2\\Ji,qi 



2J, + 1 



(30) 



Using the generalized Wigner-Eckart theorem for the 
spherical tensor operator Qa/j 



P'kmQxiiP'm'K' 



and the relation 



r^j M 



T.^'Ax,'Qx,'Prk' (31) 



PmkP'm'K' - Sjj'Skm'Pmk' 



(32) 



for projection operators one obtains for the reduced 
matrix element {Jf , qf\\Q2\\Ji, qi)- 
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{j},qf\\Q^u^q,) = jf f^'/'''iif)f;l:"'ii^) (33) 

K^Kf 

Jf 2 Ji 



where N = aj^ak is the particle number operator, and 



mq,)\NR{n)mq,)) 
($fe)|i?(l])|<i>(g,;)) 



drpvir;qj,q^;n) . (39) 



with Jf ^ 2Jf + 1, fi'^iq) ^ {-lYffiq) for < 0, 
and 



Q,^{K',K,;qf,q,) EE ($(g/)|Q2^P/A,^J$(g,)). (34) 



More details on the calculation of the reduced E2 matrix 
element are given in Appendix |Al The matrix elements 
of the charge quadrupole operator (^2^ = e J2p '"p^2p(f^p) 
are calculated in the full configuration space. There is 
no need for effective charges, and e simply corresponds 
to the bare value of the proton charge. 

Electric monopole (EO) transitions are calculated from 
the off-diagonal matrix elements of the EO operator. The 
corresponding diagonal matrix elements are directly re- 
lated to mean-square charge radii that provide signatures 
of shape changes in nuclei. The relation between EO tran- 
sitions and shape transitions and coexistence phenomena 
has been extensively investigated [1^ |35| - [37| . The EO 
transition rate t{EO) between 0'^ and can be sep- 
arated into two factors: the electronic and the nuclear 



1 



t{EO) 

where the nuclear factor pji is defined by: 
pl,m= (0+|f(i?0)|0+)'/eX 



(35) 



(36) 



and T{EO) = e^r^. The off-diagonal matrix elements 
of the EO operator can be evaluated using angular mo- 
mentum projected GCM wave functions: 

(o+|f(£;o)|o+) 

= E/o+fe)/o+(90(*(9j-)ini?0)Po"o|a>(g,)) .(37) 

Finally, it will be useful to check the average number of 
particles for a collective state \'^a'^)'- 



J2 AK,K./:''^^(g,)/i^^(gO 

qj ,qi;Ki,K2 

X J dnDilj,^mq,)\NR{n)mq,)) , 



Pvif') is the zeroth component of the nucleon vector 
current (cf. Eq. [Sc]), and the expression for the cor- 
responding transition vector density pv{r]qj,qi]Vt) has 
been given in Ref. [l^. Since the intrinsic state \^{qi)) 
corresponds to a BCS wave function, i.e. it is not an 
eigenstate of the particle number operator, the trace of 
the transition density in Eq. (j39p generally does not equal 
the total nucleon number. 



III. THE LOW-SPIN SPECTRUM OF ^^MG 

In this section we perform several illustrative configu- 
ration mixing calculations that will test our implemen- 
tation of the 3D angular momentum projection and the 
generator coordinate method. The intrinsic wave func- 
tions that are used in the configuration mixing calcula- 
tion have been obtained as solutions of the self-consistent 
relativistic mean-field equations, subject to constraint on 
the axial and triaxial mass quadrupole moments. The in- 
teraction in the particle-hole channel is determined by the 
relativistic density functional PC-Fl [2^, and a density- 
independent (S-force is used as the effective interaction 
in the particle-particle channel. Pairing correlations are 
treated in the BCS approximation. The pairing strength 
parameters Vt (t = p, n) are adjusted by fitting the av- 
erage gaps of the mean-field ground state [111 of ^^Mg 



(A) 



J2k fkvl 



(40) 



(38) 



to the experimental values obtained from odd-even mass 

differences using the five-point formula: An = 3.193 
MeV, and A^^^ = 3.123 MeV. The quantities fk arc de- 
fined in Eq. p^ and vl are the occupation probabilities of 
single-nucleon states. The resulting pairing strengths are 
V„ = 511.300 fm^ MeV for neutrons, and Vp = 518.350 
fm'^ MeV for protons. We note that these values differ 
from the universal parameters of Ref. [2^ , that have been 
adjusted to pairing properties of heavy nuclei. With the 
original pairing strengths of Ref. [2^ , the resulting gaps 
for ^''Mg are considerably smaller than the ones obtained 
from experimental odd-even mass differences. 

Parity, D2-symmetry, and time-reversal invariance are 
imposed in the mean-field calculation, and this implies 
that the space-like components of the single-nucleon four- 
currents {j^,jxv) vanish. The scalar {ps,PTs) and vec- 
tor {pviPtv) densities in the EDF of Eq. © are sym- 
metric under reflections with respect to the j/z, xz and 
xy planes. Obviously these symmetries are not fulfilled 
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by the transition densities and, therefore, the octant 
x,y,z > must be extended to the entire coordinate 
space when evaluating transition densities. 

To solve the Dirac equation for triaxially deformed 
potentials, the single-nucleon spinors are expanded in 
the basis of eigenfunctions of a three-dimensional har- 
monic oscillator (HO) in Cartesian coordinate [s^ with 
iVsh major shells. In Ref. [25| it has been shown that 
A'sh = 8 is sufficient to obtain a reasonably converged 
mean-field potential energy curve for ^''Mg. The HO 
basis is chosen isotropic, i.e. the oscillator parameters 
bx ~ by = bz = bo = ^/h/mujQ in order to keep the basis 
closed under rotations [40, |41| . The oscillator frequency 
is given by hutQ = 41A~^/'^. The Gaussian-Legendre 
quadrature is used for integrals over the Euler angles 0, 
and ip in the calculation of the norm and hamiltonian 
kernels. With the choice of the number of mesh points 
for the Euler angles in the interval [0,7r]: = N^p = 8, 
and Ng = 12, the calculation achieves an accuracy of 
« 0.05% for the energy of a projected state with angular 
momentum J < 6 in the ground-state band psl ]. 

The nucleus ^^Mg presents an illustrative test case 
for the 3DAMP-I-GCM approach to low-energy nuclear 
structure. The principal motivation for considering this 
nucleus is the direct comparison of the present analy- 
sis with the resuhs of Ref. [H, where a 3DAMP4-GCM 
model has been developed based on Skyrme triaxial 
mean-field states that are projected on particle number 
and angular momentum, and mixed by the generator co- 
ordinate method. Collective phenomena are, of course, 
much more pronounced in heavy nuclei and, therefore, 
the goal is to eventually apply the present approach to 
the rare earth nuclides and the Actinidc region. This will 
require not only a large oscillator basis, but also a large 
number of mesh-points for the Gaussian quadrature in 
coordinate space, as well as a finer mesh for the Euler 
angles and the deformation parameters. 

Axially symmetric AMP+GCM calculations are at 
present routinely performed for heavy nuclei 0] , and from 
such studies one can estimate that Nf 16 shells have 
to be included in the oscillator basis for the systems in 
the mass region around Pb. Note that the computing 
time necessary for the evaluation of one overlap matrix 
element scales approximately with Nj. For instance, 
the number of mesh-points in the axial deformation /3 
that was used in Ref. Q is a factor 4 larger than in the 
present analysis and, moreover, in the 3D case one also 
needs a finer mesh for the integration over Euler angles. 
These considerations show that a straightforward appli- 
cation of the existing 3DAMP+GCM codes to A ^ 200 
heavy nuclei will basically depend on the availability of 
large-scale general-purpose computer resources. On the 
other hand, the introduction of additional approxima- 
tions could considerably reduce the computing require- 
ments. For instance, the overlap functions are strongly 
peaked at q = q' , and the use of Gaussian overlap approx- 
imations has produced excellent results in many cases. 
These approximations form the basis for the derivation 



of a collective Bohr Hamiltonian for quadrupolc degrees 
of freedom 



A. Convergence of the 3DAMP+GCM calculations 

The convergence of the 3DAMP-I-GCM calculation has 
been examined with respect to both the number of mesh 
points in the (/?, 7) plane, and the cutoff parameter C that 
is used to remove from the GCM basis the eigenstates of 
the norm overlap kernel yy"^ with very small eigenval- 
ues n-j^/n'^^g^^ < C- In the first step the cutoff is set to 
C = 5 X 10~^, and we compare low- lying spectra of ^^Mg 
that are obtained in 3DAMP-I-GCM calculations with 
different numbers of points of the discretized generator 
coordinates. We consider the following sets of generator 
coordinates: (AI, AH, AIII) include only axial deforma- 
tions (prolate and oblate shapes) 

. AI: (/3,7) = (0.1,0°), (0.3, 0°), (0.5, 0°), (0.7, 0°), 
(0.9, 0°), (1.1, 0°); 

. AH: {13, j) ^ (0.1,0°), (0.1, 60°), (0.3, 0°), (0.3, 
60°), (0.5, 0°), (0.5, 60°), (0.7, 0°), (0.7, 60°), (0.9, 
0°), (0.9, 60°), (1.1, 0°), (1.1, 60°); 

. AIII: (/3,7) = (0,0°), (0.1, 0°), (0.1, 60°), (0.2, 0°), 
(0.2, 60°), (0.3, 0°), (0.3, 60°), (0.4, 0°), (0.4, 60°), 
(0.5, 0°), (0.5, 60°), (0.6, 0°), (0.6, 60°), (0.7, 0°), 
(0.7, 60°), (0.8, 0°), (0.8, 60°), (0.9, 0°), (0.9, 60°), 
(1.0, 0°), (1.0, 60°), (1.1, 0°), (1.1, 60°) . 

(TI, TII, Till) denote different sets with 7 7^ (triaxial 
shapes) : 

• TI: 7 = 0°,30°,60°; 

• TII: 7 = 0°, 20°, 40°, 60°; 

• TIH: 7 = 0°, 10°, 20°, 30°, 40°, 50°, 60°. 

The sets of {(3, 7) mesh-points shown in Fig. [5] have been 
used in the present analysis. 

In TableUwe display the excitation energies and B{E2) 
values of low-spin yrast states in ^^Mg, calculated with 
the 3DAMP-I-GCM model, but including only axially de- 
formed mean-field states (coordinate sets AI, AH and 
AIII, as shown in Fig. [5]). One notes that the largest dif- 
ferences in the calculated excitation energies are within 
10 % and the B{E2 : J ^ J- 2) values agree within 5 %. 
The major step for the energies comes from the inclusion 
of oblate shapes (AH) in the GCM configuration mixing 
calculations. It lowers the total ground state energy by 
« 300 keV and increases the energies by w 150 keV for 
the 2^, 200 kcV for the 4^, and 150 keV for the 6^ state. 
The refinement of the mesh in AIII produces only small 
changes. 

Similar results are found when comparing results of 
3DAMP+GCM calculations based on triaxial intrinsic 
states: AI+TI, AI+TII and AI+TIII in Table II, and 
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FIG. 2: (Color online) The distribution of mesh points in 
the {13, 7) plane for the sets AI, All, AIII, AI+TI, AI+TII, 
AI+TIII, AIII+TI, AIII+TII and AIII+TIII. 



TABLE I: Ground state energies Egs, excitation energies Ex 
(in MeV) and B{E2) values (in e^fm*) for transitions between 
low-spin states in ^"Mg, calculated with the 3DAMP+GCM 
model for the generator coordinate sets AI, All and AIII (see 
text for details). 



quantities 




AI 


All 


AIII 






-196.985 


-197.291 


-197.279 


Ex(2+) 




2.196 


2.351 


2.330 






5.394 


5.905 


5.849 


E46t) 




10.426 


10.591 


10.568 


B{E2 : 


^0+) 


78.155 


78.721 


79.135 


B{E2 : 4+ - 




137.679 


140.814 


139.750 


B{E2 : 6f - 


^4+) 


177.025 


169.246 


168.527 



AIII+TI, AIII+TII and AIII+TIII in Tabic III. The ef- 
fect of including triaxial deformations, i.e. the 7 degree 
of freedom, is perhaps best illustrated in the comparison 
between results obtained with the sets of generator co- 
ordinates AIII (Tab. U) and AIII+TIII (Tab. The 
inclusion of triaxial states in the GCM configuration mix- 
ing calculation lowers the total energies by 39 keV for 0^, 
180 keV for 2+, 226 keV for 4+, and 262 keV for 6;,^. The 
corresponding B{E2 : J — > J — 2) values are enhanced 
by 3.94%, 5.29%, 7.63% for = 2^ , if, ef , respectively. 

The influence of the 7 degree of freedom, and the con- 
vergence of 3DAMP+GCM calculations with respect to 
the number of mesh points of the discrctizcd generator 
coordinates can clearly be seen in the comparison of cal- 



TABLE II: Same as Table[Il but for the generator coordinates 
sets AI+TI, AI+TII and AI+TIII (see text for details). 



quantities 


AI+TI 


AI+TII 


AI+TIII 




-197.285 


-197.304 


-197.307 


E42+) 


2.241 


2.198 


2.177 




5.776 


5.725 


5.677 


E46+) 


10.485 


10.413 


10.360 


B{E2 : 2+ 


^ 0+) 80.523 


80.849 


81.435 


B{E2 : 4+ 


2+) 144.441 


145.926 


147.178 


B{E2 : 6+ 


-^4+) 171.275 


178.015 


182.199 


TABLE III: Same as Tabled but for the generator coordinates 
sets AIII+TI, AIII+TII and AIII+TIII (see text for details). 


quantities 


AIII+TI 


AIII+TII 


AIII+TIII 




-197.290 


-197.306 


-197.318 


E42+) 


2.239 


2.205 


2.189 


E.K) 


5.735 


5.695 


5.662 


E.{&+) 


10.452 


10.388 


10.345 


B{E2 : 2t - 


-> Oi) 80.498 


81.488 


82.256 


B{E2 : 4+ - 


2+) 143.042 


145.525 


147.137 


B(E2 : 6+ - 


4+) 166.952 


175.157 


181.379 



culations with mean-field states at the mesh points of 
coordinate sets AI, AI+TIII and AIII, AIII+TIII. The 
inclusion of triaxial shapes lowers the energies by « 300 
keV. On the other hand, very similar results are obtained 
in calculations based on coordinate sets that differ only 
in the number of axial points. Therefore, we find that, 
if prolate as well as oblate configurations are included, 
the spectroscopic properties of low-spin states in ^^Mg 
are not very sensitive to the number of axial mcshpoints. 
The inclusion of the 7 degree of freedom changes this 
situation somewhat, but not dramatically for the ground 
stated band where the admixtures with K are small. 
This is consistent with the results of the 3DAMP+GCM 
calculation with particle-number projection [2^ . based 
on the non-relativistic Skyrme density functional. It was 
shown, namely, that the number of axial states that can 
be added to the set of triaxial states is not large. Re- 
dundancies appear very quickly in the norm kernel when 
more states are added to the nonorthogonal basis, and 
this is simply a consequence of very few level crossings 
as function of deformation in ^''Mg. In Table IIVI we 
show the excitation energies and B{E2) values for low- 
lying states in ^^Mg, calculated with the 3DAMP+GCM 
model based on a set of axial mean-field states with 
13 = 0,0.1,0.2, ••• ,1.1 and 7 = 0°,60°, as functions of 
the cutoff parameter C, that defines the basis of ^^natural 
states" . Eigenstates of the norm overlap kernel with 
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eigenvalues n'l/n:^^^^. < ( are removed from the GCM 
basis {n^^^^ is the largest eigenvalue of the norm kernel 
for a given angular momentum). The excitation energies 
are not sensitive to the particular value of the cutoff pa- 
rameter provided ^ < 1 x 10~^, whereas the effect on the 
B{E2) values is < 1% for smaller values of C,. However, ( 
cannot be taken arbitrarily small, because spurious states 
are introduced in the basis for very small eigenvalues of 
the norm overlap kernel. The remaining calculations pre- 
sented in this work have been performed using the value 
C = 5 X 10-3. 



I l\M « I I I I I I I I I I I I I I I I 




FIG. 3: (Color online) Upper panel: energies and the av- 
erage axial deformations for the two lowest GCM states 
with angular momentum 0^ , 2^ , 4"'" , 6^ in ^''Mg, together 
with the mean-field (dotted) and the corresponding angular- 
momentum projected energy curves. Lower panel: squares of 
collective wave functions [^^(g)^ with 522 ~ for the cor- 
responding lowest GCM states in '^^Mg. These results are 
obtained in the axial IDAMP-I-GCM calculation. Positive 
(negative) values of the axial deformation /3 correspond to 
prolate (oblate) configurations. 



B. Axially-symmetric AMP+GCM calculation 

By restricting the set of intrinsic states to axially sym- 
metric configurations: 7 = and 7 = 180°, the com- 
plicated 3DAMP-I-GCM model is reduced to a relatively 
simple IDAMP-l-GCM calculation. For the choice of gen- 
erator coordinates /? = 0, 0.1, 0.2, • • • ,1.1; 7 = and 
7 = 180°, we have calculated the energies and the av- 
erage axial quadrupolc deformations of the two lowest 
GCM states, for each angular momentum: 0+, 2+, 4+, 
and in ^^Mg, as shown in Fig. [3] 
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FIG. 4: (Color online) Lowest energy levels of angular mo- 
mentum J'" = 0+, 2+, 4+, 6+ in ^''Mg, and reduced 
E2 transition probabilities in e^fm*, calculated with the 
3DAMP+GCM and IDAMP-f GCM models. See text for de- 
tails. 



The mean-field energy surface is somewhat soft with a 
prolate deformed minimum at /3 « 0.50,7 — 0°; ^-^d f^e 
total energy E = —192.807 MeV. This result is consistent 
with our previous calculation that used the PC-Fl energy 
density functional plus a monopole pairing force [3l , and 
with an earlier study that employed the rclativistic mean- 
field model with the NL2 effective interaction [s^ . A ro- 
tational yrast band is calculated in the prolate minimum, 
with the squares of collective wave functions (probabili- 
ties) concentrated at /3 « 0.5. 

In Fig. 3] we display the lowest energy levels of an- 
gular momentum J"" = 0+, 2+, 4+, 6+ in ^^Mg, cal- 
culated with the 3DAMP-hGGM and IDAMP-hGGM 
codes, for the sets of axially symmetric generator co- 
ordinates: (3 = 0, 0.1, 0.2, • • • ,1.1 with both prolate 
(7 = 0) and oblate states (7 = 60° and 7 = 180° 
in 3DAMP+GCM and IDAMP-l-GCM models, respec- 
tively) (columns I and III), and with only prolate 
states 7 = (columns II and IV). As expected, the 
3DAMP+GCM and IDAMP-hGCM calculations pro- 
duce virtually identical results, with small differences at- 
tributed to the numerical accuracy. In fact, the differ- 
ence between the B{E2) values shown in columns I and 
III can be further reduced by increasing the number of 
mesh-points used in the Gaussian-Legendre quadrature 
over the Euler angles (f>, 9 and ip in the calculation of the 
norm and hamiltonian kernels. 



C. Triaxial AMP+GCM calculation 

In Fig.[5]we plot the self-consistent RMF+BCS triaxial 
energy surface of ^^Mg in the /3-7 plane (0 < 7 < 60"), 
obtained by imposing constraints on the expectation val- 
ues of the quadrupolc moments (720 and (722- The panel 
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TABLE IV: Excitation energies (in MeV) and B{E2) values (in e^fm*) for transitions between low-spin states in ^'*Mg, 
calculated with the 3DAMP+GCM model for the generator coordinates 13 = 0, 0.1, ■ • ■ , 1.1, and 7 = 0°, 60° as functions of the 
cutoff parameter ( that defines the basis of '"natural states'^ . 

C 5 X 10"^ 1 X 10-2 5 X 10"^ 1 X 10-3 5 x 10"^ 1 x 10"^ 

E:,{2t) 2.275 2.340 2.330 2.341 2.314 2.320 

E^{4+) 5.703 5.931 5.849 5.573 5.544 5.580 

B{E2 : 2+ -> 0^) 77.586 80.138 79.135 79.967 79.113 79.688 

B{E2 : 4+ 2+) 144.403 143.314 139.750 136.372 137.989 138.669 
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FIG. 5: (Color onhne) Self-consistent RMF-I-BCS energy sur- 
face (left panel) of ^^Mg in the /3-7 plane (0 < 7 < 60°), and 
angular momentum projected energy surface with J'" = 0"*" 
(right panel). The contours join points on the surface with 
the same energy. The difference between neighboring contours 
is 1.0 MeV. 
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FIG. 6: (Color online) Average neutron and proton pairing 
gaps of 2*Mg in the /3-7 plane (0 < 7 < 60'^). The contours 
join points on the surface with the same pairing gap. The 
difference between neighboring contours is 0.2 MeV. 



on the right displays the projected energy surface with 
J"" 0+: 



(41) 



The contours join points with the same energy and the 
difference between neighboring contours is 1.0 MeV. The 



energy surfaces nicely illustrate the effects of including 
triaxial shapes and of the restoration of rotational sym- 
metry. The mean-field energy surfaces are found to be 
quite soft with a minimum at an axial prolate deforma- 
tion /3 0.5. When compared with the axial plot in 
Fig. [3l one realizes that the oblate minimum on the ax- 
ial projected energy curve with J'^ = 0+ is actually a 
saddle point in 7 direction. Projection shifts the mini- 
mum to a slightly triaxial shape with j3 = 0.50,7 = 20° 
and E = -197.074 MeV. The gain in energy from the 
restoration of rotational symmetry is 4.266 MeV. The 
fact that angular momentum projection leads to triaxial 
minima in the PES was already noted in 3DAMP calcu- 
lations in the eighties |48| . and very similar results have 
been obtained recently [2J] for the nucleus ^'^Mg using 
the Skyrme functional SLy4. We note, however, that the 
3DAMP+GCM model used in Ref. [Ij includes a projec- 
tion on proton and neutron numbers, that is not carried 
out in the present analysis. 

Fig. [6] displays the corresponding average neutron and 
proton pairing gaps (A), defined by Eq. PO]) . as functions 
of deformation variables /3 and 7. The gaps are relatively 
small around the minimum of the potential energy sur- 
face (PES), whereas larger values are calculated at the 
saddle points. The fluctuations of pairing gaps reflect the 
underlying shell structure. 

The solution of the HWG equation ([T4| yields the exci- 
tation energies and the collective wave functions for each 
value of the total angular momentum and parity J'^ . In 
addition to the yrast ground-state band, in deformed and 
transitional nuclei excited states are usually also assigned 
to (quasi) /3 and 7 bands. This is done according to 
the distribution of the angular momentum projection K 
quantum number (Figs. |81 IT0|) . Excited states with pre- 
dominant K = 2 components in the wave function are 
assigned to the 7-band, whereas the /3-band comprises 
states above the yrast characterized by dominant K = 
components. As an example, in Fig. [7] we display the 
low-spin PC-Fl excitation spectrum of ^^Mg obtained 
by the IDAMP-fGCM calculation with the AIII set of 
generator coordinates, and by the 3DAMP+GCM calcu- 
lation with the AIII-I-TIII set of mesh points, in compar- 
ison with available data j45l - l47| . The level scheme is in 
rather good agreement with data, but in both cases the 
calculated spectra are systematically stretched as com- 
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FIG. 7: (Color online) The low-spin level scheme of ^''Mg calculated using the 3DAMP+GCM model and IDAMP+GCM 
model with the PC-Fl relativistic density functional, in comparison with data [45l - [47| . The B{E2) values are given in units of 
e^fml 
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FIG. 8: (Color online) Contour plots of the probability dis- 
tributions Ig'af for the ground state 0^ and the first excited 
state 2^ (both the K = and K — 2 components) in ^"^Mg. 



pared to experimental bands. This is because angular- 
momentum projection is performed only after variation 
and, therefore, time-odd components and alignment ef- 
fects are neglected. Cranking calculations, for instance, 
correspond to an approximate angular-momentum pro- 
jection before variation j49| . and lead to an enhance- 
ment of the moments of inertia in better agreement with 
data [53, [5l|. However, at present the full 3D angular- 
momentum projection before variation, plus GCM con- 
figuration mixing, is still beyond the available computing 
capacities. The agreement of the calculated quadrupole 
transition probabilities with data in Fig.[7]is remarkable, 
especially considering that the calculation of B(E2) val- 
ues is parameter-free, i.e. the transitions are calculated 
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FIG. 9: (Color online) Same as in Fig. [S] but for the excited 
states 3^ and 2^ in ^*Mg. 



employing bare proton charges. 

In Fig. [51 we plot the corresponding distributions l^^p 
of Eq. with respect to /? and 7, for the ground 

state of and the first excited state 2^ (both the K = 
and K = 2 components) in ^^Mg. These quantities give 
the probabilities that the intrinsic wave functions of the 
corresponding states have a certain quadrupole deforma- 
tion characterized by the collective coordinates /3 and 7. 
For the ground state, and for the K = component of 
2^ , these distributions are largely concentrated along the 
prolate symmetry axis. Since the K = component of 
the 2^ state exhausts 92% of the norm, this state obvi- 
ously belongs to the K = band built on the nearly pro- 
late ground state. From the PES shown in the right panel 
of Fig. O with the pronounced minimum at 7 w 20°, one 
would have expected the maximum of the probability dis- 
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FIG. 10: (Color online) Same as in Fig. [3 but for the excited 
states 4^ and 4^ in ^*Mg. 
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FIG. 11: (Color online) Average particle numbers for the 
3DAMP+GCM states belonging to the three bands of ^*Mg 
in Fig. [71 



tributions in this region of the (/3, 7) plane. However, it 
turns out that the inclusion of quadrupole fluctuations 
through GCM configuration mixing, drives the structure 
built on the ground state back toward the prolate sym- 
metry axis, i.e. the GCM model calculation docs not 
predict the existence of a stable triaxial structure of the 
intrinsic states of the ground-state band of ^''Mg. The 
probability distributions for the excited states 3i and 2 J 
are shown in Fig. [HI For 2^ state the K = 2 component 
exhausts about 87% of the norm and, therefore, 2^ and 
3^ are assigned to the K = 2 (quasi) 7 band. The K = 0, 
K = 2 and K — A probability distributions of the states 
4^ and 4^ are displayed in Fig. (TU] Since the A' = 



{K = 2) component of the state 4^ (4^) exhausts 92% 
(79%) of the norm, 4'^ belongs to the ground-state band, 
and 4^ to the (quasi) 7 band. 

Finally in Fig. [TT] we display the average neutron and 
proton numbers (cf. Eq. for the 3DAMP-fGCM 

states belonging to the three bands of ^^Mg in Fig. [T] 
The dispersion of the particle number is relative large 
(w 0.3) for states 0^ and 2+. 



IV. SUMMARY AND OUTLOOK 

The framework of relativistic energy density function- 
al has been very successfully applied to the description 
of a rich variety of structure phenomena over the whole 
nuclear chart. However, to go beyond the modeling of 
bulk nuclear properties and perform detailed calculations 
of excitation spectra and transition probabilities, one 
must extend the simple single-reference (mean-field) im- 
plementation of this framework, and include long-range 
correlations related to restoration of symmetries broken 
by the static mean field and to fluctuations of collective 
coordinates around the mean-field minimum. Building 
on recent models [1, @ that have employed the genera- 
tor coordinate method (GCM) to perform configuration 
mixing of axially-symmetric relativistic mean-field wave 
functions, and especially on Ref. [l^, where we have al- 
ready considered three-dimensional angular- momentum 
projection (3DAMP) of relativistic mean-field wave func- 
tions, in this work a model has been developed that uses 
the GCM in configuration mixing calculations that in- 
volve 3D AMP wave functions, generated by constrained 
self-consistent mean-field calculations for triaxial nuclear 
shapes. 

The current implementation of the relativistic 
3DAMP+GCM model has been tested in the calcula- 
tion of spectroscopic properties of low-spin states in 
^''Mg. Starting with the relativistic density functional 
PC-Fl [2^, and a density-independent J-force as the 
effective interaction in the pairing channel, the intrin- 
sic wave functions arc generated from the self-consistent 
solutions of the constrained RMF+BCS equations in 
the basis of a three-dimensional harmonic oscillator in 
Cartesian coordinates. The constraints are on the ax- 
ial and triaxial mass quadrupole moments. After restor- 
ing rotational symmetry by 3DAMP, the fluctuations 
of quadrupole deformations are included by performing 
GCM mixing of angular-momentum projected configu- 
rations that correspond to different values of the gen- 
erator coordinates /? and 7. The GCM calculation has 
been tested both with respect to the number of mesh- 
points in the discretized (/?, 7) plane, and the cutoff- 
parameter that is used to eliminate from the GCM basis 
the "high momentum" eigenvectors of the norm overlap 
kernels with extremely small eigenvalues. Results for ex- 
citation energies in the ground-state, (quasi) 7 and /3 
bands, and the corresponding intcrband and intraband 
transition probabilities have been compared with avail- 



13 



able data on low-spin states in Mg. The comparison 
has shown a very good agreement between data and the 
predictions of the relativistic 3DAMP-I-GCM model. 

The choice of ^^Mg allows a direct comparison of 
the present analysis with the results of Ref. [2J|, where 
a 3DAMP-I-GCM model has been developed based on 
Skyrme triaxial mean-field states that are projected on 
particle number and angular momentum, and mixed by 
the generator coordinate method. Because it includes 
projection on particle number, the model of Ref. [IJ] is 
much more involved and the numerical implementation is 
more difficult. In particular, the use of general EDFs in 
GCM calculations, i.e. energy functionals with an arbi- 
trary dependence on nuclcon densities, leads to disconti- 
nuities or even divergences of the energy kernels as func- 
tions of deformation, that can possibly produce spurious 
contaminations in the calculated excitation spectra (for a 
detailed discussion, we refer the reader to Refs. j53 - l54l |. 
and references cited therein). Even though the results 
of the present calculation for ^^Mg are in good agree- 
ment with those of Ref. , an important advantage of 
performing particle-number projection is that it prevents 
a collapse of pairing when the level density around the 
Fermi energy is reduced as, for instance, close to the min- 
imum of the potential energy surface. The comparison 
with Ref. [l^l thus points to an obvious improvement of 
our 3DAMP-I-GCM model, i.e. the implementation of 
particle-number projection. 

As an alternative approach to five-dimensional 
quadrupole dynamics that includes rotational symmetry 
restoration and takes into account triaxial quadrupole 
fluctuations, one can construct a collective Bohr Hamil- 
tonian with deformation-dependent parameters. In a 
recent work (42j . we have developed a new implemen- 
tation for the solution of the eigenvalue problem of a 
five-dimensional collective Hamiltonian for quadrupole 
vibrational and rotational degrees of freedom, with pa- 
rameters determined by constrained self-consistent rela- 
tivistic mean-field calculations for triaxial shapes. As in 
the present work, in addition to the self-consistent mean- 
field potential of the PG-Fl relativistic density functional 
in the particle-hole channel, for open-shell nuclei pairing 
correlations are included in the BCS approximation. In 
Ref. H^, the model has been applied in the study of 
shape phase transitions in the region Z = 60, 62, 64 
with N « 90. The collective Hamiltonian can be derived 
in the Gaussian overlap approximation (GO A) Q to the 
full five-dimensional GCM. With the assumption that the 
GCM overlap kernels can be approximated by Gaussian 
functions, the local expansion of the kernels up to second 
order in the non-locality transforms the HWG equation 
into a second-order differential equation for the collec- 
tive Hamiltonian. Therefore, having developed both the 
five-dimensional quadrupole collective Hamiltonian, and 
the full 3DAMP-I-GCM model, we plan to perform mi- 
croscopic tests of the GOA in a study of low-spin spec- 
troscopy of 7-soft transitional nuclei, especially the effect 
of GOA on the calculated transitions between bands. In 



general, we expect that both models will be a useful ad- 
dition to the theoretical tools that can be used in studies 
of complex structure phenomena in medium-heavy and 
heavy nuclei, including exotic systems far from stability. 
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Appendix A: Reduced matrix element of the 
quadrupole operator 

The basic expressions for the calculation of EM tran- 
sition probabilities in the framework of an AMP-I-GCM 
approach are given in Ref. 0. Here we start from the 
reduced matrix element of the quadrupole operator in 
Eq. (j33p and derive a formula for the overlap matrix el- 
ements Eq. 

Q2^.iK',K;q,,qj) = ($(g,)|Q2^PiJ,^|<I>(g,)) 

= J dnD'k*,KmQ2^Rmv, (Al) 

with the overlap function of the quadrupole operator 

{Q2^.m))^J = mq^)\Q2^RmmqJ)) 

= TT[Q2^p''m{R{n)),,. (A2) 

The expressions for the norm overlap {R{il))ij and tran- 
sition densities p*-'(ri) are given in Eqs. (A28) and (C4) 
of Ref. 0. 

The indices qi,qj run over all generator coordinates. 
For Hq points on the coordinate mesh, only nqiriq + l)/2 
overlaps need to be evaluated, for instance those with 
qi ^ qj- The remaining part with qi > qj is determined 
by simply exchanging the indices i and j: 

Q2^{K',K-q^,q,) (A3) 
= ^IJ^E / dnDi^,{n)Dl,^mQl,^R{mr 

In the derivation of above relation, an irreducible tensor 
Q2-^i been introduced as = (~l)^Q2-^i- 
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The matrix elements of the multipole moment operator angular momenta: 
Qxfj, = r^Yxfj, in the spherical harmonic oscillator basis 
\nljm) read: 

(Qa^)™™' = {nl\r^\n'l') ■ {ljm\Yx,,\l'fm') . (A4) {ljm\Yx^\l' j'm') = {-1^ 
The radial part is given by 

inl\r^\n'n= (-1)"-' [r^FK)] 



J A / 



-TO fi m 



Olinii/) , 

(A6) 



E 



[r(n + Z + i)r(n' + Z' + i)]i/ 

r(t + fj) 



/ , 7, I /2 ^ ("^5) where the irreducible matrix elements of the spherical 



harmonic are given by the expression 



(cr - l)!(n - cr)!(n' - a)\{(T + n)\{a + v' - n)\ ' 



with the integers t = \{l + V + X + I), v = - l + 
and v' = \{l~V + X). 

Apart from parity conservation (Z + T + A = even), 
the angular part does not depend explicitly on orbital 



Oi|yA||/> = (-ir 



/ji'A / ] A f 



in 



1 i 

2 2 



(A7) 
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